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The development of hybrid quantum systems is central to the advancement of emerging 
quantum technologies, including quantum information science and quantum-assisted sen- 
sing. The recent demonstration of high-quality single-crystal diamond resonators has led to 
significant interest in a hybrid system consisting of nitrogen-vacancy centre spins that 
interact with the resonant phonon modes of a macroscopic mechanical resonator through 
crystal strain. However, the nitrogen-vacancy spin-strain interaction has not been well 
characterized. Here, we demonstrate dynamic, strain-mediated coupling of the mechanical 
motion of a diamond cantilever to the spin of an embedded nitrogen-vacancy centre. Via 
quantum control of the spin, we quantitatively characterize the axial and transverse strain 
sensitivities of the nitrogen-vacancy ground-state spin. The nitrogen-vacancy centre is an 
atomic scale sensor and we demonstrate spin-based strain imaging with a strain sensitivity of 
3xlO~^ strain Hz~^/^. Finally, we show how this spin-resonator system could enable 
coherent spin-phonon interactions in the quantum regime. 
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Nitrogen-vacancy (NV) centres in diamond are among the 
most promising implementations of quantum bits for 
quantum information processing^ and nanoscale field 
sensors^. Because of their excellent coherence properties^ and 
their abiHty to coherently couple to various external fields, NV 
centres are ideal candidates for integration into a hybrid quantum 
system, which combines the merits of two or more physical 
systems to mitigate their individual weaknesses. To date, 
experimental efforts have focused on integrating the NV centre 
into hybrid systems consisting of photons , nuclear spins^ and 
magnetically coupled mechanical resonators^'^. 

Recent demonstrations of high-quality (Q) factor single-crystal 
diamond (SCD) mechanical resonators^^'^^ have generated interest in 
the development of a monolithic structure consisting of such a 
resonator that is strain coupled to the electronic spin of an embedded 
NV centre^^'^^. With the NV's long quantum coherence time and 
diamond's low-mechanical losses, this system provides a unique 
environment to engineer spin-spin entanglement or spin 
squeezing^ ^. Moreover, a coherent spin-phonon interaction may be 
used to generate non-classical states of the resonator, such as 
Schrodinger cat states^^"^^. An inherent advantage of a strain- 
coupled spin-resonator system is the intrinsic nature of the coupling; 
no additional components are required to tune the coupling strength, 
such as an atomic force microscope tip or cavity, which may 
introduce fluctuating fields that decohere the NV spin. In addition, in 
contrast with other hybrid mechanical systems, this NV-resonator 
system requires no functionalization of the resonator with, for 
example, a mirror or electrode, which could adversely affect the Q. 

However, to engineer such a system, we require a quantitative 
understanding of how strain couples to the ground- state spin 
levels of the NV centre. Despite the NV's well-understood 
sensitivity to magnetic^ and electric ^^'^^ fields and 
temperature the NV ground-state strain sensitivity has not 
been fully characterized, and there exist competing theoretical 
models of the strain interaction 

In this letter, we use high-Q SCD cantilevers to probe the 
interaction between a single NV spin and crystal strain. The 
magnitude of the strain is highly controlled, and is set by the 
cantilever's mode shape and amplitude of motion, the depth and 
lateral location of the NV in the cantilever and the orientation of 
the NV with respect to the cantilever axis. The strain induced by 
an oscillating cantilever modulates the energy of the ground- state 
spin levels, which is detected via coherent quantum control of the 
spin^^. Combined with the known strain profile of a cantilever 
under pure bending and the four possible NV orientations in the 
diamond crystal structure, our measurements allow us to extract 
the axial and transverse strain susceptibilities. Thus, the NV is a 
calibrated strain sensor, with which we image the strain profile of 
the resonator. Importantly, we demonstrate NV-strain coupling 
in a geometry appropriate for reaHzing a hybrid quantum system, 
and show that by scaling down the dimensions of the resonator, it 
should be possible to reach the quantum regime of coupling 
between spins and phonons. 

Results 

Strain model and experimental setup. Here, we model strain by 
an effective electric field in which strain -induced displacements of 
the atoms alter the electron density in the crystal, resulting in 
local electric fields^^'^^. Therefore, the ground-state Hamiltonian 
for an NV in the presence of magnetic and strain fields may be 
described by (/z = 1) 

Hnv = {Do + d\\e\\)Sl + y^yS-B-^{e-'^^S\ +e'^^Sl) 

(1) 



where Dq is the crystal-field splitting (Dq = 2.87 GHz), y^y is the 
gyromagnetic ratio (7nv = 2.8MHzG~ ^), d\\ and d_i are the 
strain susceptibility parameters parallel and perpendicular to the 
NV symmetry axis (whose values are determined in this work), 
S± are the spin-1 raising and lowering operators, e\\ — e^y 

e± = + e^, tan(0J = Cy/e^ and {ei}i^x,y,z the diagonal 

components of the strain tensor defined in the NV's basis (note 
we have neglected shear). We will henceforth refer to d\\ and d± 
as the axial and transverse strain susceptibility parameters. A 
uniaxial strain along the NV symmetry axis, ey , results in a shift 
of the crystal-field splitting, whereas a strain transverse to the 
symmetry axis, e± , results in mixing between the | ms = ± 1 ) Sz 
eigenstates. Because strain and stress are tensors, the uniaxial 
stress due to the bending of the cantilever will generate strain 
along all three principal directions because of the Poisson effect^^. 
This effect cannot be neglected, as shown later. 

In this experiment, an NV interacts with the fundamental 
mechanical mode of the cantilever via strain (Fig. la). For small 
displacements of the cantilever, the strain felt by the NV centre is 
linear in cantilever position, and the total strain can be written as 
e = 60 ^, where is the amplitude of driven motion, Xq is the 
amplitude of zero point motion and eo is the strain induced by the 
zero point motion^^. If the NV is in the presence of a static 
magnetic field aligned closely to its symmetry axis, equation (1) 
can be approximately diagonalized with eigenstates {|0), | 
and the energies can be expressed as a function of the normalized 
beam displacement, X — ^ 

Eo ^ 0; £±(X) ^ Do+^||X ± (y^yB.f + (g^xf (2) 

where — d\\eo\\ and g± — d±eo± are the axial and transverse 
strain couplings to the zero point motion of the resonator, and 
I ± ) are defined as (Supplementary Note 1), 

|0) = |0); |±)=cos^|l) Te-'^^sin^l-l) (3) 

where a/2 is the Stiickelberg angle and tan(a) = d^e^/y^yB^. We 
note that by tuning we may selectively characterize axial and 
transverse strain, as evidenced in equation (2). 

The experiment consists of a single spin embedded in a SCD 
cantilever (Fig. lb). The cantilevers used here (Fig. Ic) were 
fabricated using the technique described in ref. 10, and are 60 |im 
long, 15|im wide and ~ 1 |im thick. They have fundamental 
frequencies oj^/Ik near 1 MHz and quality factors exceeding 
3.0 X 10^, which were measured using optical interferometry . 
All measurements were performed at room temperature and at 
pressures of approximately 10~^Torr. NVs were formed 
approximately 50 nm below the diamond surface via ^"^N ion 
implantation and vacuum annealing (see Methods). The 
orientation of the NV with respect to the cantilever 
(Supplementary Fig. 1) is identified through its Zeeman 
splitting in an applied static magnetic field. We align this field 
closely to the NV axis to break the degeneracy of the | ± 1) spin 
levels. We encode our qubit in the |0) and | + ) spin states, and 
coherently manipulate the spin using near- resonant microwaves. 
We are able to distinguish the effects of ey and e^by noting an AC 
ey modulates the qubit frequency at cOm, whereas an AC e± 
modulates the qubit frequency at approximately 2(jj^. 

Axial strain detection. Figure 2a shows the measurement pro- 
tocol for axial strain detection. The cantilever is resonantly driven 
with a piezoelectric transducer, and the phase of the drive is not 
phase-locked to the timing of the control pulses. The NV spin is 
initialized to the |0) state via optical pumping with a 532 nm laser. 
Microwaves tuned to the |0) ^ | + ) transition are applied to the 
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Figure 1 | A hybrid NV-cantiiever system, (a) Strain modulates the spin levels of an NV centre ennbedded in a cantilever. The strain profile was simulated 
using a finite element method, (b) Schematic of experiment: a green laser is focused onto a single NV spin embedded in a cantilever via a confocal 
microscope for initialization and readout. Microwaves are used for pulsed spin manipulation, (c) Confocal image of a cantilever showing the presence of 
single embedded NV centres. 



NV to carry out a Hahn echo pulse sequence with a total free 
evolution time T—lx. During the free evolution periods, the 
resonantly driven motion of the cantilever imprints a relative 
phase onto the qubit, which is then converted into a population 
difference and read out via the spin-dependent fluorescence. In 
the measurements of axial strain, was set to 22 G to suppress 
the effects of transverse strain, as shown in equation (2). 

A dynamic axial strain results in a modulation of the spin 
population difference that is correlated to the motion of the 
cantilever. As depicted in Fig. 2b, when the total evolution time of 
the spin is equal to an odd integer multiple of the cantilever period, 
the accumulated phase is maximized, resulting in a reduced 
population in |0). When the total evolution time is equal to an 
even multiple of the cantilever period, the qubit acquires no phase, 
resulting in a revival of |0) state population. Because the cantilever's 
motion and microwave pulse timing are not phase locked, the signal 
averages over a uniformly distributed initial phase of the motion, 
and is given by a zero order Bessel function (Supplementary Note 2). 

By fitting the data in Fig. 2b (solid black line) to the 
expected echo signal, we extract an axial strain coupling 
G|| =^I|X= 2.1 ± O.lMHz and a cantilever frequency 
cOm/271 = 882.0 it 2.0 kHz, in agreement with the driving 
frequency of the cantilever of 884.583 kHz. The strain coupling 
can be tuned with cantilever drive. Figure 2c shows a Hahn echo 
measurement of a single NV for different beam deflection 
amplitudes. In the absence of cantilever drive, the spin retains its 
coherence over a single period of the cantilever's motion. As the 
amplitude of driven motion increases, the increased strain 
coupling results in a reduced |0) state population. For sufficiently 
high couplings, Bessel fringe oscillations are observed, indicating 
that the spin has precessed more than once around the equator of 
the Bloch sphere. In Fig. 2d, we plot Gy for four different drive 
amplitudes and confirm that the coupling is linear with beam 
displacement. 

Single spin strain imaging. W^ith its atom- sized spatial extent, 
the NV is a novel, nanoscale strain sensor. W^e demonstrate NV- 
based strain imaging by measuring Gy for several NVs located at 



different positions along the length of the cantilever 
(Supplementary Fig. 2). In Fig. 3a, we show the simulated strain 
profile of a cantilever under pure bending. Figure 3b shows the 
measured strain couplings for several NVs along the cantilever for 
a fixed drive. The data are in good agreement with the theoretical 
strain profile, shown by the shaded region in Fig. 3b 
(Supplementary Note 3). Notably, this agreement provides con- 
vincing evidence that strain is responsible for the spin evolution. 
From the data in Fig. 3b, we can extract d\\ by combining the 
interferometrically measured amplitude of motion, the NV's 
position and orientation in the cantilever, and the expected strain 
profile of the cantilever (Supplementary Note 4). By measuring 
several NVs, we average over the variations in NV depth and local 
strain inhomogeneities, and we find an average value of 
= 13.4 ± 0.8 GHz per strain, and hence calibrate our strain 
sensor. Given our current experimental parameters, we estimate 
an axial strain sensitivity of approximately 3xlO~^ strain 
Hz~^^^ (Supplementary Note 5). For the cantilevers in this 
experiment, an NV at the base could then detect ~7nm of 
motion in a 1 sec measurement. 

Transverse strain detection. We next measured the transverse 
strain coupling. To do so, we measured NVs oriented perpendi- 
cular to the cantilever stress, as shown in Fig. 4a. Although the 
stress is entirely perpendicular to NVs oriented [III] and [ill] 
(shown with red bonds), these NVs will experience non-negligible 
axial strain in addition to transverse strain because of the Poisson 
effect^^. Experimentally, the effects of axial and transverse strains 
are identified by beatnotes present in the spin evolution. To 
enhance our sensitivity to these beatnotes, we applied an XY-4 
control sequence^^ (Fig. 4b) to the NV and decreased to 16 G. 
This sequence increases our sensitivity by extending the 
coherence time, correcting for first- order timing errors, and 
increasing the interrogation time. The data and fit to the expected 
XY-4 signal (Supplementary Note 6) in Fig. 4c give 
G^ =6.7 ±0.4 MHz, Gy = (2.7 ± 0.5) x 10^ kHz and a cantilever 
frequency C0m/27i = 899 ± 8 kHz, in agreement with the 
cantilever drive frequency of 884.890 kHz. Transforming the 
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Total evolution time 7"(|lis) 




Total evolution time 7"(|lis) Beam deflection (nm) 

Figure 2 | Axial strain detection with a single NV. (a) AC axial strain is detected with a Hahn echo sequence. Transverse strain is suppressed by an 
applied 22 G axial magnetic field. The n pulse allows the phases from the first (light grey) and second (dark grey) free evolution periods to add 
constructively, (b) Measured Hahn echo spin population (red circles) for an NV at the base of a cantilever with ft>m/27r = 884 kHz and an oscillation 
amplitude of 650 nm as a function of total evolution time 7= 2t. The echo signal is given by a zero-order Bessel function (fit is black solid line) because the 
cantilever motion is not phase-locked to the timing of the control pulses. Error bars correspond to standard error in photon counting, (c) The amplitude of 
cantilever motion impacts the spin evolution. Hahn echo spin population of an NV at the base of a ft»m/27r = 630 kHz cantilever in the absence of drive 
(green triangles), with a drive corresponding to 150 nm cantilever amplitude (blue squares) and 300 nm amplitude (yellow circles). Dotted black lines are 
fits to the data, (d) Plotted strain coupling G for an NV 11 [xm from the base on the cantilever described in b as a function of beam deflection. Vertical error 
bars correspond to standard error of the fit to the theoretical echo signal. Horizontal error bars correspond to the uncertainty in beam deflection from 
interferometric measurements. 



cantilever strain tensor into the NV's coordinate system and 
using the Poisson ratio of 0.069 for CVD diamond^^ 
(Supplementary Note 4), we extract ^:/|| = 13.3 ± 1.1 GHz 
per strain and = 21.5 ± 1.2 GHz per strain. The extracted 
value of is in very good agreement with the extracted value of 
d\\ measured for the NVs experiencing predominantly axial strain 
in Fig. 3. 

Discussion 

The measured strain susceptibiHty parameters indicate that while 
the simple model equating strain and electric field accurately 
predicts the qualitative nature of the NV-strain coupling, it fails 
to capture the quantitative nature of the strain susceptibility 
parameters. In particular, this model predicts a value of d\\ that 
is approximately 50 times smaller than dj^ (ref. 18), whereas 
our measurements show that d\\ is approximately 1.5 times 
smaller than d±. During the preparation of this manuscript, 
a closely related study^^ reported measurements of 
^^11 =5.46 ±0.31 GHz per strain and = 19.63 ± 0.40 GHz per 
strain using static bending of diamond cantilevers. In comparison 
with our results, we find good agreement with their value of d± 
and find that their value of d\\ is approximately three times 
smaller than our measured value. We make several observations 
in regards to this discrepancy. First, Teissier et al?^ report large 
variations (up to one order of magnitude) in the strain 
susceptibility parameters across NVs. Second, the AC sensing 
techniques used in our work are inherently more sensitive than 
the techniques used in ref. 27, and hence, the amplitude of 
cantilever bending is significantly smaller here. Large bending 
amplitudes could lead to deviations from Euler-BernouUi theory. 



Furthermore, the resonance phenomenon we utilize provides a weU- 
defined mode shape, whereas static bending may not isolate a single 
mode. Another study^^ measured shifts of the NV ground-state 
crystal field splitting under hydrostatic pressure, and their measured 
value of d\\ — 17.5 GHz per strain is in good agreement with our 
result. However, no value oi d^ was reported. Ah initio calculations 
from ref 23 estimate a value of (iy = 7.44 GHz per strain by 
considering modifications of the spin-spin coupling because of 
changes in the NV bond lengths. This value is of the same order of 
magnitude as our result, indicating that this effect has a significant 
role in the NV-strain interaction. 

Looking towards quantum appHcations^^, our cantilever's 
current zero point motion coupling of — 0.04 Hz may be 
significantly increased by miniaturizing the resonator. For 
instance, a doubly clamped beam with dimensions of (2 jim, 
50 nm, 50 nm) and 09^/271 = 238 MHz has a coupling 
^ 180 Hz. Although having smaller strain couplings than, for 
example, quantum dot systems where the coupling can be 
hundreds of kHz (ref. 28), the NV centre benefits from its long T2 
coherence times. Assuming T=100mK, a mechanical Q=10^ 
and 2L T2— 100 ms, this doubly clamped beam achieves a single 
spin cooperativity, rj — ln^^ ^ 1.6, where y — cOm^Q is the 
damping rate and n is the average thermal occupation number of 
the mechanical mode. We note that although smaller resonators 
would produce a larger gj_, clamping losses would limit their 
Qs^^. At T<50K, SCD resonators have demonstrated^^ Q> 10^. 
In the devices described here, the Hahn echo T2 times are roughly 
10 |is. Using NVs formed by a nitrogen delta-doping technique^^ 
and applying higher order dynamical decoupling, longer T2 times 
are possible and recently, a T2 > 500 ms has been measured at 
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low temperatures . A cooperativity f/ > 1 enables experiments in 
the quantum regime that are mediated by a coherent spin- 
phonon interaction^^. We note that the strain coupling to the 
orbital levels of the excited state is significantly higher"*'^^, and for 
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Figure 3 | Single-spin strain imaging, (a) Stress profile for cantilever used 
in this experiment at 250 nm of beam deflection using a finite element 
method simulation, (b) Measured strain coupling (orange circles) as a 
function of the NV's distance from the cantilever base for a fixed 250 nm 
oscillation amplitude. The grey shaded area shows the region of expected 
strain couplings from theory including uncertainties in NV depth (13 nm) and 
amplitude of driven motion (10 nm). Vertical error bars correspond to the 
standard error from fits to the expected control sequence signal. Horizontal 
error bars correspond to the uncertainty in the NV's lateral position. 



the resonator described above, would give a coupling gj_ of 
several MHz, which presents an alternative route to accessing the 
quantum regime. 

In conclusion, we have presented a novel hybrid system in 
which a single NV centre spin is strain -coupled to the mechanical 
motion of a macroscopic cantilever. We quantitatively measured 
the NV ground- state spin coupling to axial and transverse strain, 
and these results will guide theoretical investigations of the strain 
coupling mechanism. Although our measurements indicate that 
shear does not couple significantly to the ground- state spin levels 
of the NV centre, we note that the cantilever-based technique 
presented here can be extended in future experiments to probe 
effects of shear on the NV. Furthermore, the coherent strain 
detection approach we have presented can be extended to other 
solid-state spin systems, such as SiC colour centres^^'^"*. With its 
atomic-scale spatial resolution, NV-based strain sensing could be 
useful for identif)^ing the poorly understood origins of dissipation 
in mechanical resonators, such as anchoring losses, which will be 
instrumental in engineering higher Q resonators^^. Moreover, our 
results will provide insight into the role of strain- induced 
decoherence that could ultimately limit the effectiveness of 
many NV-based applications^^'^ Most importantly, we 
demonstrated that our system provides a feasible route to 
accessing coherent spin-phonon interactions, with applications 
ranging from quantum control of a macroscopic resonator to 
long-range spin- spin interactions. 

Methods 

Experimental setup and sample preparation. The experimental setup consists of 
a home-built confocal microscope with an integrated vacuum chamber. The 
experiments are done at room temperature and high vacuum (~ 10~^Torr) to 
reduce mechanical dissipation in the resonators. A continuous wave laser at 
532 nm is used for optical pumping and readout of the NV spin, and is gated with 
an acousto-optic modulator (AOM). Photons emitted by the NV into the phonon 
sideband are collected into a single-mode fibre and directed to a fibre-coupled 
avalanche photodiode. For strain measurements, the degeneracy of the ms = ± 1 
spin levels is broken by a DC magnetic field supplied by a rare earth magnet. 
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Figure 4 | Measurement of transverse strain, (a) Stress profile for a cantilever during bending (red arrows). The face of the cantilever (gold square) 
corresponds to the (110) plane. NVs oriented along [111] and [111] (red bonds) experience mostly transverse strain, with a snnall axial strain because of the 
Poisson effect^^. NVs oriented along [111] and [111] (green bonds) experience predonninantly axial strain (see Supplementary Information). Measurements 
of transverse strain are done with NVs oriented [111] and [111], (b) XY-4 pulse sequence used to measure transverse strain. Transverse strain (in blue) 
modulates the qubit splitting approximately twice as fast as axial strain (in red), (c) XY-4 spin population (see Supplementary Information) for an NV 
(yellow circles) 3 \irr\ from the base with a beam deflection of 675 nm. The expected XY-4 signal for a model considering both axial and transverse strain is 
shown in black. Error bars correspond to photon counting noise. 
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Microwaves used for resonant spin manipulation are delivered to the sample by a 
wire that sits approximately 30 \xm above the sample. Phase control of the 
microwave pulses is accomplished via IQ modulation. Timing of the pulse 
sequences is controlled by a Spincore Pulseblaster ESR-Pro 500 MHz card. 

The cantilevers used in this experiment were fabricated using the techniques 
described in ref 10. NVs were formed by ^'^N ion implantation with a dosage of 
3 X 10^ ion cm ~ at 40keV and a 0° tilt, which yields an expected depth of 
51.5 nm (calculated by Stopping and Range of Ions in Matter (SRIM)). The sample 
was annealed under high vacuum ( ~ 10 ~ ^ Torr) at 800 °C for 3 h. 

The cantilevers are resonantly driven with a piezoelectric transducer that is 
clamped to the sample holder and electrically driven with a lock-in amplifier 
(Zurich Instruments). The cantilevers' mechanical properties are characterized 
using standard optical interferometry techniques, as described in ref 10. For our 
experiments, it is very important that we know the amplitude of driven motion of 
the cantilever. This is calibrated by finding the piezo voltage at which the tip of the 
cantilever hits the silicon substrate that sits 1 |im below the cantilever. This is 
measured by an interruption of the Bessel fringes in the interferometric signal. 

Error analysis. Errors in the extracted axial strain coupling parameters and 
resonator frequencies are given by the standard error in the fit where the couplings 
and the resonator frequency are the only free parameters. The error in the quoted 
axial strain susceptibility parameter is given by the standard error in the average 
value of those parameters as measured by nine NVs. The error in the transverse 
strain coupling and susceptibility parameter are given by the measurement of a 
single NV, and the standard error of the fit shown in Fig. 4c. 
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